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Transitions in turbulent rotating Rayleigh-Benard convection 
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Numerical simulations of rotating Rayleigh-Benard convection are presented for both no slip and free slip boundaries. The goal is to 
find a criterion distinguishing convective flows dominated by the Coriolis force from those nearly unaffected by rotation. If one uses heat 
transport as an indicator of which regime the flow is in, one finds that the transition between the flow regimes always occurs at the 
same value of a certain combination of Reynolds, Prandtl and Ekman numbers for both boundary conditions. If on the other hand one 
uses the helicity of the velocity field to identify flows nearly independent of rotation, one finds the transition at a different location in 
parameter space. 
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^ 1 Introduction 
Oh 



At least two flow regimes exist in rotating Rayleigh-Benard convection in plane layers: Near the onset 
of convection, the nonlinear advection terms are small and the Coriolis force dominates the dynamics, 
provided the Ekman nmnber is small enough. The controlling effect of the Coriolis force is the distinguishing 
feature of the first regime. If the Rayleigh number is increased at fixed Ekman and Prandtl numbers, the 
advection term in the Navier-Stokes equation becomes larger than the Coriolis term (which is linear 
in velocity) so that rotation becomes irrelevant and the flow behaves as if rotation was absent. This 
deflnes the second regime . For the purpose of predicting the Nusselt number, it was found useful in 
Schmitz and Tilgner to introduce a transitional regime into the classification, in which the Nusselt 

number obeys a power law different from those observed in the two other regimes. There is an ongoing 
debate concern i ng the p a ramet e rs at which th e tran si tion betwee n the fi r st and second regime s occu rs 
(|Rossbvl (|l969l ). lAurnoul (|2007l ). iLiu and Eckl (|2009l ). iKing et all (|2009l ). ISchmitz and Tilgnej (|2009l )). 
While there is by now ample evidence that the naive criterion that the Rossby number equals one at the 
tra nsition is inadequ ate, there is no agreement on what the correct criterion is. Recently, it was shown 
by King et~al\ ( 20091 ) that experimental heat flux data are compatible with the idea that the flow is in 
one regime or the other depending on wh ether the thermal bo u ndary layer is thicker than the Ekman 
layer or vice versa. A numerical study by ISchmitz and Tilgne 3 (|2009l ) avoided Ekman boundary layers 
by employing stress free boundary conditions. A transition between the two regimes still occurs and an 
analysis of the asymptotic behavior of the Nusselt number leads to a transition criterion based upon a 
combination of the Reynolds, Prandtl and Ekman numbers. 

The present paper follows up on the numerical simulations presented in ISchmitz and Tilgner (|2009l ) 
and inte nds to answer two question s: First, how important are the boundary conditions? Is the criterion 
found in Schmitz and Tilgner 

m) stiU relevant for no slip boundaries? And second, is the classification 
of a given flow into t he di ffere nt regimes indepe ndent of the quantity used for that classification? Both 
Schmitz and Tilgner (I2OO9I ) and lKing et aH (|2009l ) classified fiows according to their Nusselt number. Here, 
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we will also look at a quantity derived from the velocity field which is important for the dynamo effect, 
namely the helicity. 

2 The mathematical model 

A plane layer of thickness d, filled with fluid of kinematic viscosity u, thermal diffusivity k, and thermal 
expansion coefficient a rotates with angular velocity Q about the z— axis. This axis is perpendicular to the 
layer. Gravitational acceleration g is pointing in the negative z— direction. The temperatures of the top 
and bottom boundaries are fixed at Tq and Tq + AT, respectively. These two boundaries are assumed to be 
either free slip or no slip, whereas periodic boundary conditions are applied in the x— and y— directions. 
The equations of evolution are made non-dimensional by using cP/k, d and AT for units of time, length, 
and temperature, respectively. These equations then become within the Boussinesq approximation for the 
dimensionless velocity v{r,t) and temperature T(r,t): 

Pr 

dtv + (v ■ V)v + 2 z xv = -Vp + Pr V^v + Ra Pr Tz (1) 

Ek 



V ■v = (2) 



dtT + v-S/T = V^T (3) 

z is the unit vector in z— direction and p collects the pressure and the centrifugal acceleration. The 
boundary conditions state in terms of the adimensional quantities that T{z = 0) = 1, T{z = 1) = 0, 
and free slip boundary conditions require that Vz = dzVx = dzVy = at both z = and z = 1, whereas 
Vx = Vy = Vz = for no slip boundaries. Three independent dimensionless control parameters appear: The 
Rayleigh number Ra, the Ekman number Ek, and the Prandtl number Pr. They are defined by: 

Ra=- , Ek = —^ , Pr = -. 4 

The Reynolds number Re and the Nusselt number Nu are computed as: 



The overbar denotes average over time and the integrals extend over the computational volume V for Re 
and over the surface A of either the top or the bottom boundary for Nu 



The equations of motion were solved with the same spectral method as used in ISchmitz and Tilgner 



Resolutions reached up to 129 Chebychev polynomials for the discretization of the z— coordinate 



and 256 x 256 Fourier modes in the (x, y)— plane. The periodicity lengths along the x— and y— directions 
were always chosen to be identical. Since the typical size of flow structures varies considerably as a function 
of the control parameters in rotating convection, it is not useful to use a single aspect ratio throughout 
all simulations, where the aspect ratio is defined as the ratio of the periodicity length in the plane 
and the layer height. The aspect ratio was adjusted for each Ek to fit at least 8 columnar vortices along 
both the X— and ?/— directions at the onset of convection, and kept constant as Pr and Ra were varied. 
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Figure 3. (Nu — l)Ek^^^ as a function of Re Pr Ek^^^ for the same data and with the same symbols as in figs. [Tl and l3l The dashed 
lines are power laws with exponents 2 and 2/3. The points inside the interval delineated by the two vertical lines reappear in figure [J] 
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Figure 4. Nut, as a function of Raf^. The symbols have the same meaning as in fig.[T] This figure contains only those data points which 
lie in the interval marked by vertical lines in fig.|3] The straight lines indicate the functions Nut = 0.11 ■ RoP^'^^ (lower, dot dashed line) 
which fits well to simulations with free slip boundaries, and Nu^ = 0.17 ■ Ra^/^^ (upper, dashed line). 



3 No slip boundaries 



Computations for no slip boundaries are much more demanding than in the free slip case. That's 
why the parameter range covered by the data in figure [T] is smaller than in the equivalent figure in 
Schmitz and Tikned (|2009l ^. Nonetheless, the data suffice to show that the behavior is essentially the same 



for both types of boundary conditions. Figure [U shows Nu{Ra) for no slip boundary conditions. One rec- 
ognizes the well known pattern of a delayed onset of convection followed by a ste ep rise. Nu in the ro tating 
flow may even exceed its value for zero rotation for no slip boundary conditions ( Zhong et al. ( 20091 )). The 
delayed onset of convection is also visible in the graph of Reynolds number versus Rayleigh number (figure 

El). 

Figure [3] us es a data reduction rnotivat ed by the results for free slip boundary conditions . Some of the 
arguments of Schmitz and Tilgneil ( 20091 ) are recapitulated here to make the paper self-contained. Low 
values of Re Pr Ek}^'^ correspond to laminar flows near the onset of convection. Forming the dot product 
of eq. ([1]) and v, integrating over the whole volume and averaging over time, one finds 



e = {Nu - l)Ra 



(6) 



where (- = y j < {diVj){diVj) > dV is the adimensional average dissipation rate of kinetic energy. In a 
laminar flow, one expects e oc {Re Pr)'^ / , where A is a characteristic length scale of the flow. For Pr > 
0.676, convection starts at a cri tical Rayleigh number Rac obeying Ruc oc Ek~^^'^ and forms stationary 
cells of size Ac with Ac oc Ek (IChandradekhad (|l96ll l). Eq. ® becomes Ek~'^/^Re^Pr'^ oc {Nu - l)Ra. 
Close to onset, Ra ~ Rac and therefore {Nu — 1) oc Re^Pr^Ek^^^ . For large values of Re Pr Ek^^'^ , one 
should approach the non-rotating case and find a law independent of Ek. There is no reliable theory 
for convection far from the onset so that we have to rely on numerical results. Simulations show that 
{Nu — 1) oc {Re Pr)'^^^ for zero rotation. Both asymptotes, near onset and far from it, become straight 



lines in a logarithmic plot of {Nu — l)Ek 



1/3 



vs, 



Re Pr Ek^/"^. 



Figure Oshows {Nu-l)Ek^ as a function of Re Pr Ek^^"^ for no slip boundary conditions and should be 
compared with figure 2 of lSchmitz and Tilgner The dashed lines have the slopes of the asymptotes 

which {Nu — \)Ek^/^ has to follow either in the limit of large Rayleigh numbers (when rotation plays no 
role) or near the onset of convection. The prefact ors in these power laws are o btained from best fits to the 
simulations with free slip boundary conditions of lSchmitz and Tilgner (|2009l ) so that a direct comparison 
is possible with the results of computations with no slip boundary conditions shown by points in the figure. 
It is seen that both boundary conditions have asymptotes with the same exponents. Most importantly, 
for both boundary conditions, the crossing of the two asymptotes near Re Pr Ek^l"^ = 2 separates the 
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convection dominated by rotation from the convection virtually unaffected by rotation. 

It is noteworthy that in the no slip case of figure O the asymptotes provide us wit h reasonably good fits 
throu ghout the entire parameter range. In the case of free slip boundary conditions in Schmitz and Tilgneil 
the interval 1/2 < Re Pr Ek}^"^ < 10 (indicated by vertical bars in figure [3]) required special 
treatment. It turned out to be useful to introd uce nondimensiona l parameters iii dependent of molecular 
diffusivities already discussed in previous work (IChristensenI (jlooS) , E urnoul (120071 )): Nu* = NuEk/Pr = 
Q / {pCpATQd) and Raj^ = Nu^RaEk'^/Pr = (gaQ) / {pCpQ^ cP) , where p denotes density and Cp heat 
capacity. An envelope to the data in the interval 1/2 < Re Pr Ek^/"^ < 10 is given by Nu^ = 0.11 • Ra^ff 
for free slip boundaries. Figured] shows that a similar conclusion holds for no slip boundaries after a change 
in the prefactor: Nu^ = 0.17 • Ra^jf^ is more appropriate in this case. 



4 Free slip boundaries: Helicity 

In order to quantify the influence of rotation on the flow structure, we investigate the helicity H, or more 
precisely the correlation between vorticity and velocity, defined by: 

H = (7) 

|V X v\ \v\ ^ 

This correlation is now averaged over horizontal planes, and ideally over time. However, helicity was not 
saved during the production runs so that it had to be deduced from single snapshots. Assuming ergodicity 
of the flow, the spatial average replaces to some extent the temporal average. The noise in the figures 
below is low enough so that further simulations did not seem warranted. The average helicity {H), defined 
by 

{H) = jfdxfdyH (8) 



where A is the area of a horizontal plane, is shown in figure [5j As expected ( Chandradekhar ( 196 ll )). 



rotation introduces helicity of one sign in one half of layer and of the opposite sign in the other half 
{z = 0.5 is not an exact plane of symmetry in figure [5] because of the missing time average). As Ra 
increases, the advection term increases compared with the Coriolis term and the helicity decreases. The 
maximum of {H) is reached at a point at a smaller and smaller distance from the boundaries as Ra is 
increased. Helicity is induced by the presence of the boundaries because they force in- and outflow out of 
departing or arriving plumes which is spun up by the Coriolis force. When turbulence tends to destroy the 
correlation between vorticity and velocity present in helical structures, it is easiest to do so away from the 
walls. 

A suitable global measure for the influence of the rotation on the velocity fleld is the rms fluctuation in 
{H){z), i.e. 



Hrms= [\Hy'dz. (9) 
^0 

Figure [6] gives Hrms as a function of Re Pr Ek^/"^ . There is a reasonable collapse of the data for a given 
Pr, which shows that the combination Re Pr Ek^^^ captures the right Re and Ek dependences, but the 
points for different Pr do not fall on top of each other. However, if one plots Hrms as a function of Re Ek^^"^ , 
one again obtains a good collapse of the data points, at least at large ReEk^/"^ (see figure [7]). Within the 
noise on Hrms, the fiow is indistinguishable from nonrotating convection for ReEk^^'^ > 50. The transition 
criteria are thus not identical when based on the Nusselt number or on the helicity. The transition occurs 
at a certain value of Re Pr Ek^/"^ in as far as heat flux is concerned, and at a certain value of Re Ek^/"^ for 
helicity. 
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Figure 5. Average helicity (H) as a function of height z for no slip boundaries, Ek = 2 X 10 and Pr = 7. In going from the curve 
drawn with a continuous line of largest amplitude to the one with smallest amplitude, the Rayleigh number is Ra = 1.0 X 10^, 2.0 X 10^, 
5.0 X 10•^ 1.0 X 10'', and 1.0 x lO'^. The dashed curve is for Ra = 1.0 X 10*. 
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Fi gurG 6. Hfuis 3^ ^ function of FIg. Pv Ek^^^ . The filled symbols are for Pv — T and the empty symbols for Pv — O.T. The Ekman 
numbers are Ek = 2.0 X 10"^ (diamonds), 2.0 X 10"^ (squares), 2.0 X 10"" (triangles) and 2.0 X lO"'^ (stars). The remaining symbols 
indicate no slip boundaries and Ek = 1.7 X 10"^ (x), Ek = 8.7 X 10"-* (+), Ek = 4.7 X 10'"* (*), and Ek = 2 X 10'^ (•), all for Pr = 7. 
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Fi gure T. H-rms cis a function of RcEk^^^. The symbols are the same as in figure 161 
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5 Discussion and conclusion 

It was shown in section [3] that according to the Nusselt number, Re Pr Ek}^'^ = 2 separates flows dominated 
by the Coriohs force from flows unaffected by rotation for both types of boundary conditions, free shp and 
no shp. The combination Re Pr Ek^^^ can be recognized as the Peclet number based on the thickness of 
the Ekman layer. This interpretation is not directly useful, however. It is not clear why this Peclet number 
should matter for the heat flux, because the heat flux is crossing diffusively the velocity boundary layer in 
which velocity and temperature gradient are mostly perpendicular to each other so that v ■ VT = and v 
does not enter the temperature equation ([3]) within the boundary layers. In addition, there is no Ekman 
layer for free slip boundaries, and yet, the transition occurs at the same values of the control parameters. 

On e may argue that a length scaling as Ek^/"^ is not completely foreign to free slip boundaries, either. 
Hidd (Il964l ) calculates a velocity boundary layer thickness in Ek ' near a free slip boundary provided 



that the density of the fluid varies along the surface. This cannot be the case within the Boussinesq 
approximati on on an isotherma l boundary, so that this mechanism cannot create a length in Ek^/"^ in our 
simulations. Julien et al. 1 (|l996 l) also find an Ek^/"^ layer in the velocity field if the thermal boundary layer 



thickness varies laterally. Their calculation posits a lateral temperature variation which is not advected 
by the horizontal velocity, so that this result is not directly applicable to self-consistent simulations. 
Nonetheless, these analytical calculations are motivation enough to scan the velocity profiles obtained 
from the numerical simulations for a length scaling as rapidly as Ek^^"^ as a function of the Ekman 
number. However, no such length could be found in the case of free slip boundaries. 

As section S] has shown, the transition from a flow dominated by rotation to a flow unaffected by rotation 
is not well defined, anyway. The transition occurs at different values of the control parameters depending 
on whether the classification is based on the Nusselt number or on the helicity of the velocity field. The 
helicity is down to negligible mag nitude for ReEk^'"^ > 50. We wih now apply this result to the Earth's 
core. Most data in this paper were obtained for free slip boundaries, but the close agreement between free 
slip and no slip boundaries found here encourages us to apply the transition criteria to the Earth's core 
nonetheless. For the Earth's core, the generally accepted values of 1/ = 5 X 10-^mVs, a = 7.29 x IQ-^s'^ 
and a typical fiow velocity of 5 x lead to Re Ek^/"^ = 82. This estimate casts doubts on the picture 

of a geodynamo operating as an a^— dynamo with an a— effect due to the helicity of the fiow. 
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